Meson- Baryon Form Factors in Chiral Colour Dielectric Model 

S. C. Phatak 



Institute of Physics Bhubaneswar 751 005 



The renormalised form factors for pseudoscalar meson-baryon coupling are computed in 
chiral colour dielectric model. This has been done by rearranging the Lippmann-Schwinger 
series for the meson baryon scattering matrix so that it can be expressed as a baryon pole 
term with renormalized form factors and baryon masses and the rest of the terms which 
arise from the crossed diagrams. Thus we are able to obtain an integral equation for the 
renormalized meson-baryon form factors in terms of the bare form factors as well as an 
expression for the meson self energy. This integral equation is solved and renormalized 
meson baryon form factors and renormalized baryon masses are computed. The parameters 
of the model are adjusted to obtain a best fit to the physical baryon masses. The calculations 
show that the renormalized form factors are energy-dependent and differ from the bare form 
factors primarily at momentum transfers smaller than 1 GeV. At nucleon mass, the change 
in the form factors is about f 0% at zero momentum transfer. The computed form factors 
are soft with the equivalent monopole cut-off mass of about 500 MeV. The renormalized 
coupling constants are obtained by comparing the chiral colour dielectric model interaction 
' Hamiltonian with the standard form of meson-nucleon interaction Hamiltonian. The ratio 

of ANii and NNn coupling constants is found to be about 2.15. This value is very close to 
the experimental value. 
PACS numbers: 12.39.Ki, 12.39.Fe, 13.75.Gx 
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I. INTRODUCTION 

The meson-baryon form factors are interesting and useful quantities for several reasons. For one thing, 
we expect the form factors to have some bearing on the underlying structure of the hadrons. For example, 
in perturbative chiral quark models |Q the pion-baryon form factors are directly related to the quark wave- 
functions in baryons. Form factors are essential in the effective models of meson-baryon interactions since 
hadrons, after all, are not point particles. One also needs these form factors for a consistant description of 
nuclear phenomena. The meson exchange nucleon-nucleon potentials include phenomenological vertex form 
factors which are presumably related to the structure of hadrons. Form factors are also needed in compu- 
tation of photo- and electro- production of baryon resonances since, in principle, the photon can couple to 
the virtual charged meson cloud in a baryon. Another technical reason for introducing the form factors is 
that the effective meson-baryon interaction models are, generally, not renormalizable and the form factors 
provide the needed cut-off functions. 

Usually, the form factors employed in the meson-baryon interaction models are of monopole type ( ~ 
i-fc^/Aii ) with k 2 being the four momentum carried by the meson. It must be noted that the monopole form 
is, as such, purely phenomenological and is possibly not related to the underlying structure of the hadrons. 
The form factors used in earlier calculations Q were 'hard' with the cut-off parameters A > lGeV. There 
arc, however, recent calculations H which indicate that softer form factors (A ~ 500 MeV) are required 
to fit A production on nuclei. Some indication of the value of cut-off parameter A is available from low 
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energy phenomenology. For example, Thomas and Holinde j| argue that the 3% descripancy between ppir 
and pnn + coupling constants [^) is essentially due to the four-momentum variation of the ttNN form factor 
between q 2 = m 2 and q 2 = 0. With monopole parametrization of the form factor, A of about 800 MeV is 
required to explain this descripancy. Similar arguments have been proposed by Coon and Scadron ||| for 
A ~ 800Me^. One should however note that these analyses are restricted to relatively small values of \q 2 \ 
( ~ (14QMeV) 2 ) and therefore cannot determine the form factors at larger momentum transfers. Recently, 
Saito and Afnan have investigated the dependence of triton binding energy on the nNN form factor. 
They They compute the three-body ir — it force contribution to the triton binding energy. They claim that 
one can determine the ttNN form factor their model and it turns out that their renormalized form factor is 
softer with A ~ AOOMeV. On the other hand, Schultz and Holinde ^ fit irN scattering data and obtain a 
form factor having A ~ 800AT eV. Some attempts have been made to determine the pion-baryon form factor 
from QCD. For example, Liu et al. || have used quenched lattice QCD calculations and extracted A of 750 



MeV. T. Meissner 10 , on the other hand, has used QCD sum rules and obtained A of about 800 MeV. 

In the present work, the renormalised pseudoscalar meson baryon form factors are computed by using 
the chiral colour dielectric (CCD) model of baryons O]. The CCD model is based on the calculation of 
Nielsen and Patkos Jl2| who showed that 'coarse graining' of QCD the Lagrangian on lattice gives rise to 
an effective scalar, colour-neutral field called colour dielectric field. In the CCD model, the the interaction 
of the colour dielectric field with the quark and gluon fields is such that it gives rise to the confinement 
of these objects. The chiral symmetry is restored in the CCD model by introducing the interaction of the 
pseudoscalar meson octet with the quarks [pL~T|] . The CCD model has been used to compute properties of 
light baryon [jll| as well as ttN scattering phase shifts |L3| . Generally the agreement between the CCD model 
results and the experimental data is good. In particular, a good agreement with the nN scattering phase 
shifts is obtained in the CCD model. In the ttN work, the authors showed that solving of the relativistic 
Lippmann-Schwinger equation for pion-nucleon scattering leads to the renormalization of baryon masses 
and the pion-baryon veritces are also renormalized due to multiple scattering. However the renormalized 
vertex functions ( or form factors ) were not computed in that work. The present work is essentially an 
extension of the earlier work but here we have concentrated on renormalized form factors and have not 
compute ttN scattering matrix. The basic idea here is as follows. The CCD model ( or, for that matter 
any other quark-based model with mesons coupling to quarks ) give bare meson-baryon form factors and 
bare baryon masses. These are then dressed by meson-baryon interactions. As we shall show later in the 
paper, computation of meson-baryon scattering ( including the pole positions of stable particles ) leads to 
the determination of dressed baryon masses as well as dressed meson-baryon form factors. This is essentially 
achieved by rearranging the terms in the Lippmann-Schwinger series for the meson-baryon scattering matrix 
so that one obtains equations for the dressed baryon masses ( positions of the poles in the scattering matrix 
) and meson-baryon form factors ( see the following sections for details ). Thus, one can compute the dressed 
pion baryon form factors. One should note here that these dressed ( or renormalized ) form form factors, 
and not the bare form factors should be used in nuclear physics calculations. Some saliant features, of the 
results of our calculation are 

1. The renormalization of the form factors depends on the energy of the meson-baryon system. At zero 
momentum transfer, the ratio of renormalized and bare form factors varies from about 1.1 to about 
1.3 as the energy is increased from nucleon mass to A mass. 

2. The renormalized form factors differ from the bare ones for meson momenta of 1 GeV or smaller. There 
is practically no change at larger meson momenta. 

3. The form factors computed in the CCD model are soft. If one wants to fit the CCD model form factors 
by a monopole form in the region of 300-700 Mev, corresponding A is about 500 MeV. 
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At this stage we would like to emphasize that our prescription of computing the renormalized form factors 
should be applicable to models having basic baryon-meson interaction as an input. That is, given a model 
for baryon-meson interaction Lagrangian with bare baryon masses and form factors, one should compute 
the renormalized form factors by fitting the data ( cross sections, phase shifts, positions of the poles of 
the scattering matrix etc). These renormalized form factors should then be used in other nuclear physics 
calculations such as NN potentials, photo- and electro- production of mesons etc. In other words, the 
parameters of the basic baryon-meson model are necessarily the bare parameters and should not be obtained 
directly from the experiments. If one does this then one is essentially choosing dressed parameters in the 
Lagrangian and then one should use the model at tree level and should not compute meson loops for meson- 
nucleon scattering etc since that would amount to double counting. 

The paper is organised as follows. In Section 2, a brief description of the CCD model is given and the bare 
masses and pion-baryon form factors are computed. In Section 3, begining with the Lippmann-Schwinger 
series, the the integral equations for the meson-baryon form factors are obtained. In Section 4, the results 
of the calculation are presented and discussed. 

II. CHIRAL COLOUR DIELECTRIC MODEL 

The CCD model is described and the bare masses of baryons and pseudoscalar meson-baryon vertices are 
computed in this section. The description of the CCD model is somewhat brief because in this work we want 
to focus on the computation of meson-baryon form factors. For the details of the CCD model the reader is 
refered to 

A. Lagrangian 

The Lagrangian density of CCD model is given by pd[ | 

C(x) = Mx){i$- (mo + ^rC/ 5 ) + 9 -Kf (x)}^ + ^TridpUVU*) - 

-\ml<f> 2 {x) - \x\x){F^{x)f + \al(d,x{x)? U( X ) (1) 

where U = e ,A »^ ^ and f/ 5 = e lX ^ 75 /^, ip(x), A^x), and X (x) an(1 0(#) are quark, gluon, scalar ( colour 
dielectric )and meson fields respectively, m and are the masses of quark and meson, / is the pion decay 
constant, F ilv (x) is the usual colour electromagnetic field tensor, g is the colour coupling constant and A° and 
are the usual Gell-Mann matrices acting in colour and flavour space respactively. The flavour symmetry 
breaking is incorporated in the Lagrangian through the quark mass term (mo + ^U^), where mo = for u 
and d quarks. So masses of u, d and s quarks are m, m and mo + m respectively. The meson matrix consists 
of a singlet 77, triplet of ir and quadruplet of K. 

The self interaction U (x) of the scalar field is assumed to be of the form 

U( X ) = aB X 2 {x)[l-2{l-2/a)x{x) + (l-Z/a)x 2 {x)\ (2) 

so that U(x) has an absolute minimum at x — an d a secondary inimum at % = 1. 
We would like to make the following observations on the Lagrangian of eq(l). 

1. The CCD model Lagrangian is invaraint under chiral transformations if the mass terms of the pseu- 
doscalar mesons are dropped. 
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2. The coupling of the scalar field (%) to the quark and gluon fields is such that the quark mass becomes 
infinite and the dielectric function ( x 4 ) vanishes when \ becomes zero. This means that the quarks 
and gluons cannot exist in a region where x becomes zero. That is the confinement of quarks and 
gluons is included in the model. 

3. The two minima of the self-interaction of the scalar field (eq(2)) can be identified with the perturbative 
and physical vacua of the MIT bag model Thus, the CCD model is a dynamical model generating 
the bag. 

We follow the cloudy bag model approach jjj in the present work and treat the gluon and meson interactions 
with the quarks perturbatively. Therefore we expand the Lagrangian in powers of 1// and keep the terms 
up to order 1// in the Lagrangian of eq(l). With this approximation the CCD model Lagrangian becomes, 

C(x) = 4>(x){i$- (m + -^t(i + - f \a^ a {x) + |a q 4»U 
I X{x) i 2 J 

+ \(d,Mx)) 2 ~ \™%<?{x) - \x\x){F^)? + \°l{d»x{x)f - U( X ) (3) 

The parameters of the CCD model are quark masses ( m and mo ), the 'bag constant' ( B ), strong 
coupling constant ( a s — g 2 /47r ), pion decay constant ( / ), glueball ( or dielectric field ) mass ( mcB ) and 
the constant a in U(x)- Of these parameters, value of a is chosen to be 24 since from our earlier calculations 
fllf the results are not sensitive to it. To bigin with, we choose the experimental value of the pion decay 
constant ( / = 93M eV ) in our calculations. As we shall see later, the value of / required to fit pion-nucleon 
coupling constant is very close to 93M eV. Thus we are left with five free parameters to be adjusted. In our 
earlier calculations [jllj we found that a reasonably good fit to the baryon masses is obtained for m and B 1 ^ 4 
ranging between 100 and 140 MeV. We therefore choose m and B 1 / 4 in this range and adjust itlcb, Q s and 
mo to fit nucleon, A and A masses. Computed masses of other octet and decuplet baryons are within few 
percent of their experimental masses. 



B. Bare Baryon States 

The equations of motion of quark, gluon and colour dielectric fields can be obtained from the Lagrangian 
of eq(3). Since we shall be treating the gluon interactins perturbatively, the quark and colour dielectric field 
equations of motion are solved self-consistantly. The gluon as well as pseudoscalar meson interactions are 
then treated perturbatively. This approach is similar to the one followed in the cloudy bag model calculations 
jij. Thus the equations of motion for the quark and dielectric field are, 

(cj - m - -^-r)g l {r) = -fi(r) - -/;(r) 
X(r) r 

{e i -m -^)f i {r)=g' i {T) (4) 

and 

X"(r) + -X'{r) - ^x(r) [l - 3(1 - 2/a) X (r) + 2(1 - 3/«) X 2 (r)j + £ -^-(gf(r) - ff{r)) = 0. (5) 

We have assumed spherical symmetry in obtaining these equations. The equations of motion of quark and 
dielectric fields are solved self-consistantly with the boundary conditions that as r — * oo, x( r )> ffi( r )> fi( r ) ~ * 
0, x( r ) an d 9i( r ) ~~ * constant and fi(r) — > as r — > 0. The bare baryon states are constructed by putting 
three quarks of appropriate flavour in lsi/2 orbital. Thus the baryon wavefunction is a product of symmetric 
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space wavefunction, symmetric spin- flavour wavefunction and antisymmetric colour wavefunction [O. The 
mass of this state is computed by evaluating the matrix clement of the Hamiltonian 

H Q = f d 3 x[Y,^(-^-V + -+mo)y + ^((Vx) 2 + Tl 2 ) + U(x)\ (6) 
J f X 

in the bare baryon state. Here \& is an annihilation operator of a quark in the state computed in eq(4) and II 
is the momentum conjugate to the dielectric field x- Method of coherent states [jl5| has been used to better 
account for the energy associated with the dielectric field. To this the colour-magnetic interaction between 
the quarks is added perturbatively ( see ref for details). The bare mass of the baryon is then 

M B =< B(6)\H \B(6) > +E M (7) 

where |-B(0) > is the bare baryon state having zero momentum. This state is constructed by using Peierls- 
Yoccoz projection technique Jl7]|l£|]. It is convinient to express M B as 

M° B =m GB {C? +a s Ci) (8) 

where C-f depends on the number of strange quarks in the baryon and Cf depends on the spin-flavour wave 
function of the baryon. ( Cf and Gf are, of course functions of other parameters of the model. The mass 
M B defined above is the bare mass of baryon B since it does not include the meson self energy corrections. 
In order to obtain the physical mass of the baryon one should determine the position of the pole of the 
T- matrix in the appropriate meson-baryon channel fl3| . Clearly, this mass is different from M B . In our 
calculations we treat m,Q B , a s and mo/mcB as free parameters and obtain a best fit to the baryon masses. 



C. Bare Meson Baryon Vertex 

The meson baryon interaction Hamiltonian obtained from the Lagrangian of eq(3) is 

Htm = ~ / (fx-^-^ix^X^ix^iix) (9) 

In order to compute the bare meson-baryon form factors we quantize the meson fields and evaluate the 
matrix element of Hint between bare baryon states defined in the previous section. With this, the interaction 
Hamiltonian can be written as 

H mt = -. V \B' >< B\ f f k f X , < B'Mx) l5 X^(x)^ ai (k)\B > +h.c] (10) 



It is convinient to use angular momentum formalism for evaluating the spin- flavour matrix elements. The 
details of the calculation are described in Appendix A. The end result of the calculation is 

H int = -fJ2\ B ' >< B \^= < TB,Tt,tB,li\T B >,tB> >< S B , 1, 

J VIGtt^ 

(-iy +i * J -==( a ^(k) aKtyulwWk!,-* (11) 

where the quantities in angular brackets < • • • > are the usual SU(2) Clebsch-gordon coefficients, S B and 
T B are spin and isospin of barion B, s B and t B are their projections, /z is the isospin projection of meson <f> 
and a.BB'<t> are the reduced matrix elements. The phase factor i$ is defined in the Appendix A. olbb 1 ^ are 
defined in Table 1. The bare form factor u BBI( p(k) is evaluated in the 'brick- wall' frame so that the baryons 
\B > and \B' > have momenta \k/2 > and | — fc/2 > respectively. The momentum states of the baryons are 
constructed using Peierls-Yoccoz projection technique ll7|Jf^] . 
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TABLE I. The reduced matrix elements a BB i a, are given in this table. The matrix elements for K and r\ mesons 
are given in curly brackets and square brackets respectively. 





N A E 3 A E* H* Q 


N 
A 
E 

A 

E* 

n 


5[1] {-VT8} {-v/f} 78 {-7f} 

{3} [-2] -2 {1} 2 {-2} 

{-1} 7l2 v/f [2] {-5} {-75} 7f[2] {-2} 

{72} {^f} -l[-3] {-^ 2 [2] {-4} 

732 {-V 7 ?} MVS] {-v 7 ?} 

{78} 724 { A /±f [ -V8 1 | -7§ I VTO a/¥ 1 720 

{4} {-V 7 ?} {-78} [-78] {v 7 ?} V 7 ?!-^] {^0} 

{-4} {710} [-720] 



III. MESON-BARYON SCATTERING 

Given the meson-baryon interaction Hamiltonian of the previous section, one can compute the tree-level 
meson-baryon scattering matrix. It involves pole and crossed terms shown in Fig. 1. These two terms 
constitute the driving term or potential for meson-baryon scattering. To obtain a unitary scattering matrix 
and to include multiple scattering effects the relativistic Lippmann-Schwinger equation 



(12) 



is solved. In this equation the meson-baryon channel ( including the spins and flavours of the particles ) is 
represented by the subscripts and E~ ( (p) = \/P 2 + ni% + W </> T (p) • The potential V above is computed from 



the diagrams of Fig. 1. ( see Phatak, Lu and Landau |13| for details ). 
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(a) (b) 

FIG. 1. s channel (a) and u channel (b) pole terms ( direct and crossed terms respectively ) contributing to the 
meson-baryon potential 
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We would like to make several points regarding the T-matrix equation given above (eq(12)). 

1. The T-matrix given above describes the scattering of one of the octet of the pseudoscalar mesons 
from the octet or decuplet of baryons. Thus it represents resonance production, production of strange 
mesons and baryons etc. 

2. The channel energies E 1 in the Lippmann-Schwinger equation above are defined in terms of the physical 
baryon masses since these are the physically observable states. However the masses of baryons in the 
intermediate state of the pole and crossed terms are the bare masses defined in the previous section. 

3. The poles of the T-matrix in appropriate spin-flavour channel of the meson-baryon system should 
coincide with the physical masses of the baryons. This means that the parameters of the CCD model 
( basically the glueball mass rriQ B , strong coupling constant a s and mo ) are to be adjusted to obtain 
the best fit to the physical baryon masses defined by the pole positions of the T-matrix. 

Given the potential Vp a (k' , fc), one can solve the Lippmann-Schwinger equation and obtain the phase 
shifts and cross sections for meson-baryon scattering. For this purpose, it is best to expand the T-matrix 
equation in spin-isospin channels of the meson-baryon system [fl3f and solve partial wave Lippmann-Schwinger 
equation. As mention earlier, when one solves the Lippmann-Schwinger equation renormalization of baryon 
masses as well as the form factors is already included. However one cannot compute the renormalized form 
factors using this procedure. Since we are interested in computing the renormalized form factors in the 
present work, we do not want to solve the Lippmann-Schwinger equation as such. Instead, we express the 
Lippmann-Schwinger equation as a multiple scattering series. From this series we obtain integral equations 
for the renormalized form factors by summing up parts of the series. We also obtain the expressions for 
the renormalized masses. The details of this exercise are outlined in Appendix B. The result is an integral 
equation for the renormalized form factor UBB'4>(k) 

n x-^ t - OLBB"'d>'aB"'B"AOL ,-t' 

UBB^(k) = u° BB ,Jk) + V {-IYSb-Sb-Tb-Tb - 7 - TT2 S-SjL W (t b Tb"'Tb"T B ';T^) 

V5^B'» Un^PaBB^ 

r k' i dk'u° BB , t , (j ,,(k')u B , t , B , r<l> (k)u B>IBI ^(k') 
W(S B S B ,„S B „S B ,;ll) J _ ^- _ ^ (jfe/) _ ^ {mE _ — - — — (13) 



where a = ^2a+l, E° B (k) = y/k 2 + (M B ) 2 , E B (k) = y/k 2 + M%, u^(fc) = ^Jk 2 + m 2 , M B and M% are 
renormalized and bare masses of baryon B respectively, / = A$ + A$i — S B — S B > — Tb — Tb>, Sb and Tb are 
the spin and isospin of baryon B respectively and W(- ■ -)'s are the usual Rakah coefficients. The equation 
for the renormalised propagator of a baryon B is 



G(E) = — K + — fr 'E(E)G(E) 

v ' E — M B E — Ma 



E - M% - E(E) 

where the self energy H(E) is 



(14) 



WBn ^ (-l) T °'+ S e>- T °- S °-*S B >f B> f ^dku% Bl ^k)u B , B -{k) 

= ^ -J^^'^'^l ^(k)(M B -E B ,-^k)) (15) 

where A is 0, -1/2, 1/2, for tt, K, K and r\ respectively. The renormalized propagator has a pole at 
E = M B + Thus the renormalized mass of baryon M B is given by 

M B = M% + E(M fl ). (16) 
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Note that eqs(13-16) are coupled equations since renormalized masses and form factors appear in these 
equations. Thus these equations need to be solved self-consistantly. 

We can now use the renormalized form factors and propagators in the expression for the scattering matrix. 
We then get, 



where 



T Bm '(k,k';E) =T|1 B , , +T°Ak,k') (17) 



T B<t,,B><t>> k ; E ) = a BB»-$ aB "B'<t> e-m » ^ ' 



and 



?Bt,B>* (k, k') = VE^,b'4>' & k') + I d3 P- 

R" rh" J 



vE^A^p)Tg,,,, B ,Ap,k') 



, i B"<f>"\ l% iyJ- L B"4)",B' 

E - E B (p) - w+(p) 



(19) 



Here V Bl p B ,^, (k, k 1 ) is the crossed term shown in Fig 1. A diagramatic representation of the T-matrix is 
shown in Fig. 2. 
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= • • + — 
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FIG. 2. The graphical depiction of T-matrix equation. The first line corresponds to eq() of text. The second line 
depicts the integral equation for T . The thick (thin) lines represent baryons with physical ( bare ) masses and the 
broken lines represent mesons. Filled circles represent renormalized form factor. 

The interpretation of eq(17) is obvious. The first term represents the pole or resonance contribution which 
is given in terms of physical masses of resonances and renormalized form factors or vertex functions. The 
rest of the contribution is lumped into the second term which can be considered as a background. Clearly, if 
our model is extended to include higher order terms or exchanges of heavier mesons, these will predominantly 
affect the background term. 

Usually the T-matrix is computed by solving the Lippmann-Schwinger equation using a suitable numerical 
technique ( eg matrix inversion method ]h| ) or by summing the Lippmann Schwinger series numerically to 
a desired accuracy. Alternatively one can now compute the renormalized form factors and propagators and 
use them to evaluate the T-matrix. The computation of the pole term of the scattering matrix is very simple 
in this case. For the crossed term T B ,^ B ,i,(k, k') however one has to solve an integral equation (eq(19)). Our 
computations show that the crossed term of the T-matrix can be computed iteratively and the convergence 
is rapid. We also find that the dominant contribution to the scattering matrix comes from the pole term. 
One should note that the computation of the renormalized form factor already includes a large number 
of terms involving the crossed diagram of the meson-baryon potential. It is therefore not surprising that 
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Tg^ (k, k') appearing in the T-matrix equation (eq(19)) is somewhat less important, particularly near 
the resonance. 

Note that pole term of the T-matrix Tj\ B'<p'(k, k'; E) defined in eq(18) has a singularity at E = Mb" ■ 
If this singularity occures at physical energy, it means that the particle B" can decay and that should be 
reflected in the mass of the particle. That is the physical mass defined above should have an imaginary part 
corresponding to the width of the baryon B" . 

The renormalised form factors are computed by solving eq(13) above iteratively. We find that the iterative 
procedure does converge and the convergence is reached within fifteen iterations. Note that the integrand of 
the form factor integral equation has singularities due to the energy denominators and the integral should 
be regularised properly. We have used principal value prescription for this purpose. The results of the 
calculation are discussed in the following Section. 



IV. RESULTS AND DISCUSSIONS 

In this Section we shall discuss the results of our calculation. The precedure adopted for the calculation 
is as outlined below. The coupled quark and dielectric field equations are first solved and bare baryon 
states are constructed. The bare baryon masses are computed by evaluating the matrix element of the 
Hamiltonian (eq(6)) in the momentum-projected bare baryon states and including the colour magnetic 
interaction. The baryon masses are essentially functions of rriQB, a s and mo ( m and B 1 / 4 are kept fixed for 
a given parameter set ). Also, using the momentum-projected baryon baryon wavefunctions the bare form 
factors ( vertex functions ) for baryon-meson coupling are computed. These are then used in eq(13) and 



the renormalised form factors are computed. In principle, one could use matrix inversion method |19| for 
this purpose. However we find that the dimensionality of the matrices is prohibitively large because we are 
considering coupling to all the octet of pseudoscalar mesons. We therefore solve eq(13) iteratively. We find 
that the convergence is reached after 10-15 iterations. Using the renormalised form factors the renormalised 
masses are computed by using eq(16). The parameters of CCD model ( ttigb, ot s and mo ) are then adjusted 
to fit the renormalised masses to the physical masses. 

For the comparison with the experimental data we consider the ttNN and 7riVA coupling constants at 
nucleon and delta masses respectively. Conventionally, in terms of the nonrelativistic ttNN interaction 
Hamiltonian, the nNN coupling constant J^nn is defined as jjj 

H mt = i^E™ / -fit=v(k)[a ■ kr-a(k) + /i.e.] (20) 

where a and r are nucleon spin and isospin operators respectively, v(k) is the form factor defined to be unity 
when k — im^ and a{k) is the annihilation operator for pion field. Experimental value of f^NN is 0.279 ( 
fwNN = 0.078 ). Comparing the interaction Hamiltonian of eq(ll), we have f^NN = awW7r "^-^ fc ~°) 7 " ,r 
where the renormalized form factor is calculated at E = TOjvj the nucleon mass. Similarly, f^NA = 

Computations have been done for a number of parameter sets of CCD model. The results are quite 
similar for all these sets so long as ttigb is between 800 and 1200 MeV. Therefore we have chosen one of the 
representative parameter set ( tttqb — 978.6, a s — 0.438, B 1 / 4 — 122.3, m — 122.3 and mo = 210.1 ) for a 
detailed discussion. The pion decay constant / has been fixed at its physical value ( 93 MeV ). The bare and 
renormalized irNN and 7riVA form factors evaluated at nucleon delta masses are displayed in Fig 3. It is 
clear from this figure that the renormalizes form factors differ from the bare form factor at three momentum 
transfers smaller than 1 GeV. An inspection of eq(13) shows why this happens. As the momentum transfer 
increases, the first term in the denominator of integrand also increases because of increase in ci)^(fc). As 
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a result, the contribution of the integral to the form factor decreases with the increase in the momentum 
transfer. Even at momentum transfers smaller than 1 GeV, the change in the form factor is about 10 % for 
unnk and 30% for u^a-k- Thus one could conclude that the effects of renormalization are relatively minor. 



t 1 r 



S 




Momentum ( MeV ) 

FIG. 3. Renormalized and bare form factors vs momentum in units of glueball mass. The different curves are for 
unntt ( continuous line ), u% Njr ( dotted line ) and unAw ( dashed line ). u° NAw is same as u° NNvr . 

We can extract equivalent monopole form factors from the computed renormalized form factors. If, for 
example, unntt is expanded in power series near k = 0, one obtains A ~ 800Mey. This value of A is 
about same at the one obtained by various authors P|o],|8|-|Ic| . However one must note that as function of 
momentum transfer, the monopole form factor falls off more slowly than the CCD model form factor. We 
find that a better fit to the CCD model form factor is obtained by using dipole or Gausian form. In fact, 
we find that if we replace the bare CCD model form factor by a dipole or Gaussian form factor, we get 
practically the same results. If, on the other hand, one wants that the monopole form factor should be close 
to u NN ^(k) for k ~ 300 - 500MeV, A - 500MeF. 

The computed value of the irNN coupling constant is 0.254. This value is about 9% smaller than the 
experimental Jnntt- This, by itself, can be considered as a success of the CCD model. The /nnk as well as 
the baryon can be fitted by decreasing the pion decay constant / to 86 MeV ( a 7% decrease ). It may be 
mentioned here that this change in / does not affect the computed properties of baryons very much. 

Let us now consider the ratio of ttA^A and irNN coupling constants evaluated at nucleon and delta masses 
respectively. This ratio is 6\/2/5 = 1.697 for bare coupling constants, as dictated by SU(3) flavour symmetry. 
For renormalised coupling constants this ratio is about 2.15, which is very close to the experimental value of 
2.14 (^0|. We believe that this is a very important result emerging from our calculation. This result implies 



10 



that, to begin with, the quark model value of the ratio, which is about 25% off the experimental value, is 
itself a good starting point. In other words, the bare coupling constants in the CCD model themselves are in 
qualitative agreement with the experimental values. The renormalization of these quantities, which should 
be done in any case, is able to produce a very good agreement with the experimental result. Incidently, we 
would like to note that when / = 86MeV the irNN coupling constant agrees with the experimental value 
and the ratio of nNA and ttNN coupling constants still remains close to 2.15. 



1.6 




1.0 

1000.0 1200.0 1400.0 



Energy ( MeV ) 



FIG. 4. The ratio of renormalized and bare form factor at zero momentum transfer as a function of energy. The 
continuous curve is for NNir and the dashed curve is for NAtt. 

The energy dependence of f^NN and f^NA is shown in Fig. 4. The figure shows that the coupling 
constants increase with the energy in the energy range shown. 



V. CONCLUSIONS 



The chiral colour dielectric model has been used in this work to compute renormalized meson-baryon form 
factors. This has been done by writing the Lippmann-Schwinger equation in a series form and deriving the 
integral equation for the renormalized form factor. We find that the renormalization effects are primarily 
restricted to the meson momenta smaller than 1 GeV. The increase in the form factor at zero meson momen- 
tum is about 10% at nucleon mass and about 30% at A mass. The renormalized NNir coupling constant is 
about 10% smaller than the experimental value if /, the pion decay constant is chosen to be 93 MeV. The 
ratio of N An and 7V7V7T coupling constants is found to be 2.15, which is close to the experimental value. This 
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ratio for bare coupling constants is 1.697. Thus the renormalization effects yield a good agreement between 
experiment and theory. The computed form factors are soft with the equivalent monopole parameter A 
of about 800 MeV if fitting is done at small momentum transfers. However, the CCD model form factors 
decrease faster than the monopole form factor for large momenta and the CCD form factors are closer to 
Gaussian or dipole forms. 

Strictly speaking, the computations done in this work are not fully relativistic although relativistic ex- 
pressions for energies of the baryons are used. For a completely relativistic calculation, one will need to use 
relativistic propagators for baryons and use a relativistic integral equation, such as Blankenbeckler-Sugar 
equation instead of relativistic Scrodinger equation that is used in the present work. Indeed it is technically 
possible to do this, with some increase in the complexity of the problem. However we feel that the approach 
followed in the present work is reasonable since we have restricted ourselves to kinetic energies small in 
comparison with the baryon masses. Extension of the present calculation to higher energies would require 
the above-mentioned relativistic generalizations as well as inclusion of heavier baryon resonances. 

We would like to thank Prof. R. H. Landau for useful discussions. 



APPENDIX A: BARE FORM FACTORS 



The interaction Hamiltonian is 

Hint = — 



/ d 3 a;^-*(x)75A l *(x) ( /) 4 (x) 
/ J X{x) 



(Al) 



where Aj are the SU(3) GellMann matrices acting on the flavour coordinate of the quarks. In terms of the 
pseudoscalar octet fields 0,, we define the 7r, K and r\ fields as 



= <Pl, T l 
= 01,0 
- ^1/2,1/2 



T^(0i ±i( t>2) 



03 
1 

i 



(04 + «05) 



»K+ - ^1/2,-1/2- 



= —=(04 - i<t>b) 



(A2) 



'1/2,-1/2- ^ 



b K o = 0i/2,i /2 



= ^(06 -«07) 



<P V = 00,0 = <P8- 

Extra negative sign for n~ field is introduced to conform with the usual definition of vectors in spherical 
basis and is different from usual definition of pion fields. With these definitions of meson fields, 



A a 0a = (-1)^01,-^ + ^1/2^01/2,-^ + U l/2,A<01/2,- M + ^800,0- 



Quantizing the meson fields we have 



d 3 k 



87T 3 J yj2Ui{k) 



-,ik-x 



ai {k) + e- ik - s a\{k) 



(A3) 



(A4) 



where a, is the annihilation operator of meson i and a\ is the creation operator of meson conjugate to i The 
commutation relations for the meson field operators are 



[oi, M (fc),oi = (-irs^s(k-k'), 



(k),atfl v (k')] = S id 6^ v 6(k - k') 
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and [ao,o(fe)) Qo o(^')] = 8(k — k'). Computation of the matrix element of the interaction Hamiltonian between 
bare baryon states gives 



Hi, 



= E \ B ' 



>< B\ 



B,B',i 



d 3 k 



Oi(fc) -al(k)]u° BBIi {k) < B'^alB > ■% 



(A5) 



The details of computing the bare form factor u BBIi (k) is given elsewhere jljj ]. Computing the spin- flavour 
matrix element < B'\\i<j\B > and expressing the result in angular momentum formalism, we get 



Hint — 



E 



ia BB > 



d 3 k 



4/^ 3/2 J ^Mk) 



a\-v{k)-a{ (k) ki,-„u BBI<j ,(k) 



< T B ,X,tB,li\T B ',t B ' >< S B ,l,aB,v\S B ',SB' > 



(A6) 



Here the summation over spin and isospins of baryons and mesons is implied. Also (j) implies the meson type ( 
7r, K, K and rj ) and their total isospins and projections along z-axis are represented by A and fi respectively. 
The baryon spin ( isospin ) and its projection is represented by S B ( T B ) and s B (t B ) respectively. The 
reduced matrix clement a BB ',i, and the phase factor /,/, are determined by exlpicit calculation of one matrix 
element, as is usually done. The phase factor i<p is /i, 0, 1/2 — fx and for n, K, K and r/ respectively. 



APPENDIX B: RENORMALIZATION 

One can derive the equations for renormalized form factors and masses by writing the T-matrix equation 
in a Lippmann-Schwinger series and summing parts of the series. This procedure, however, is somewhat 
cumbersome. On the other hand, same conclusions can be arrived at by considering the graphical repre- 
sentation of the series. We have therefore shown the diagrams corresponding to the individual terms of the 
Lippmann-Schwinger series for the T-matrix of eq() in Fig(5). The figure shows the diagrams upto l// 6 . It 
should be clear that one can easily generalize the conclusions drawn from the analysis of these diagrams to 
the full Lippmann-Schwinger series. 



j - i . i — ' — 



— / \ — — / \ ^ ^ 

' \ / \ ' \ ' \ / \ / /\s 

1 i i, 1 I 1 I . i, 1 ( 1 t ^_ 



— ( - " ) ' v - I 1 i t ~~ - 1 (_ 



FIG. 5. Diagrams ( upto order l// 6 ) contributing to the T-matrix. Thick lines represent physical baryons, thin 
lines represent bare baryons and broken lines represent mesons 

The diagrams in hg(7) are organized as follows. The diagrams in the first line renormalize the right-most 
vertex. The diagrams in the second line renormalize left-most vertex and the diagrams in the third line 
renormalize the mass or the propagator of the intermediate baryon. The diagrams in the fourth line are the 
mixed diagrams representing renormalization of vertices as well as mass. Finally, the diagrams in the last 
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line are the rest of the diagrams which do not give rise to renormalization. Generalizing the diagrams in 
the first line to all orders, we obtain an integral equation for the renormalization of vertex ( form factor ). 
Diagramatically, this is represented in Fig(6) where the filled circle represents the renormalized form factor. 



/ / ' ~ ^ s 

/ / / v 

— — < — + — — 1 — 

FIG. 6. Diagramatic representation of vertex renormalization 

The corresponding integral equation obtained after some angular momentum algebra is 

u BB , 4) (k)=u BB , 4> (k)+ V (-iyS B „S B ,„T B „T B — — - B - B ^W(T B T B ,nT B nT B ,;T <j) T 4> ,) 

4,>^b»> 12n f aBB '* 

/k dk t^BB'" ^'{k ^^B*f' (jyik^u g /f gf-^ f {k ) 
- E%,„ ^Qp) ^(k))(E E B „ ^(fc')) (B1) 

where a = y/2a + T, E° B (k) = ^Jk 2 + (M B ) 2 , E B (k) = ^/k 2 + M 2 , uj^k) = ^Jk 2 +m 2 , M B and M% are 
renormalized and bare masses of baryon B respectively, I = A0 + A^> — S B — S B i — T B — T B i with A^ being 
1, 0, 1, for 7r, K, K and r\ respectively, S B and T B are the spin and isospin of baryon B respectively and 
W{- ■ -)'s are the usual Rakah coefficients. 

_ = _ + ^ + r " r ' + r^ 7" VN + • • • 



+- 1 * + • • 1— • + 



= — + — 

FIG. 7. Diagramatic representation of renormalized propagator 

Now consider the mass renormalization. The relevent diagrams are shown in Fig(7). Clearly, the third 
line in this diagram is the Schwinger-Dyson equation for the propagator with the self energy being given by 
the meson loop. The thing to notice, however, is that one of the meson-baryon form factor appearing in the 
self energy is a renormalized one where as the other is a bare one. Assuming that the computation is done 
in the frame in which the baryon is stationary, the renormalized and bare propagators are E _ M v_^ E ^ an d 
E ^ M o respectively. The self energy S(-E) is 

x-- (-1)T b >+s b/ -t b -s b -a s B ,f B , r fc 4 dfcu^ B , (fc)u i3 , i j^(fc) 

= S¥ -JM ""**"*"* I ^(k)(M B -E B ,-^(k)) (B2) 

where A is 0, -1/2, 1/2, for ir, K, K and r\ respectively. The mass of the physical baryon M B = M B +a(E = 
M B ). 

It should be evident that using the renormalized form factors and propagators the T-matrix can be written 
in the form given in eqs(16-18) and Fig(2). 
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